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Abstract 



We calculate near zero recoil the order a s corrections to the Bjorken and 
Voloshin sum rules that bound the B — * Z)M £[> form factors. These bounds 
are derived by relating the result of inserting a complete set of physical states 
in a time ordered product of weak currents to the operator product expansion. 
The sum rules sum over physical states with excitation energies less than a 
scale A. We find that the corrections to the Bjorken bound are moderate, 
while the Voloshin bound receives sizable corrections enhanced by A/Aqcd- 
With some assumptions, we find that the slope parameter for the form factor 
h^x in B —* D*£u decay satisfies 0.4 < P% < 1.3. 
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I. INTRODUCTION 



It is possible to perform model independent extractions of some of the CKM matrix 
elements and quark masses from exclusive and inclusive semileptonic B meson decays via 
a systematic expansion in inverse powers of the heavy bottom and charm quark masses. 
The form factors in B — > D^*' I v decays are related by heavy quark symmetry ]I[ to the 
Isgur-Wise function, = v ■ v f ), where v is the four- velocity of the B and v' is that 
of the D^*\ A model independent determination of \V c b\ from the differential decay rate 
dT(B — > D* lv)l dw is made possible by the fact that £ is equal to unity at zero recoil (w = 1) 
IjH-Q. Inclusive B decay rates can be calculated by performing an operator product expansion 
for the time ordered product of two weak currents ||, allowing for a model independent 
determination of \V c b\ from the inclusive semileptonic decay rate, T(B — > X c iv). 

The major theoretical uncertainties in the determination of \V c b\ from inclusive decays are 
due to the questionable convergence of the perturbative corrections to the b quark decay rate 
H, and the uncertainties in the b and c quark masses. Uncertainties in the determination of 
\V c b\ from B — > D* tv originate from order AQ CD /m^ fe corrections at zero recoil, and from 
extrapolating the form factors measured at w > 1 to w = 1 (phase space vanishes at w = 1). 
The uncertainties in this extrapolation would be reduced if the slope of the form factor at 
zero recoil were known. 

Sum rules have been derived that relate the exclusive decay form factors to the inclusive 
decay rates. The Bjorken sum rule |7]|§ gives the bound p 2 > 1/4, where p 2 is minus the 
slope of the Isgur-Wise function. Voloshin derived the upper bound p 2 < 1/4 + [mu — 
mQ)/[2(mM 1 — tum)) — 0.75 @, where tum — rriQ = A is the mass difference between 
the ground state heavy meson and the heavy quark that it contains (up to corrections of 
order Aq CD /m^ b ), and — rnu is the mass of the first excited meson state above the 
pseudoscalar- vector doublet. 

The most recent experimental data from CLEO [|K| is p 2 B _ D * = 0.84 ± 0.12 ± 0.08. 
This might violate the above upper bound when experimental uncertainties decrease. The 
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ALEPH flTJ result p\ _^ D * — 0.29 ± 0.18 ± 0.12 is significantly smaller, close to the above 
lower bound. The slope of the form factor 1%a 1 , which occurs in B — > D* tv decay, has also 
been studied by CLEO |12|]. Central values for its slope parameter p\ ranging between 
0.91 and 1.53 have been obtained. Thus, it is interesting to calculate the corrections to the 
Bjorken and Voloshin bounds. 



II. REVIEW OF SUM RULES 



To derive the sum rules, we follow Refs. ||13|-[15|]. Consider the time-ordered product 



T ^ = 2^ fd'xe«~(B\T{Jl(x),J v {0)}\B), (1) 
where Ja is a b —+ c axial or vector current, the B states are at rest, q is fixed, and 



% = m B — E M — e. Here E M = \Jm? M + \q | 2 is the minimal possible energy of the hadronic 
final stateQ that can be created by the current J M with fixed \q\. With this definition of 
e in terms of the hadronic variables, the cut of in the complex e plane corresponding 
to physical states with a charm quark lies along < e < +oo. It will be important that 
at the same value of |<f| the cut at the parton level lies within the smaller region e > 
A(w — l)/w + O (Aq CD / m c,b) ■ (T^ u has another cut corresponding to physical states with 
two b quarks and a c quark that lies between —2Em > e > — oo. This cut will not be 
important for our discussion.) To separate out specific hadronic form factors, one contracts 
the currents in (p]) with a suitably chosen four-vector a, yielding 

a- T M o" = ^-E W ft) W/WW-'W + . . . , (2) 

2m B x £j X — £j M — e 

where the ellipses denote the contribution from the cut corresponding to two b quarks and 

a c quark. The sum over X includes the usual phase space factors, i.e., d 3 p/2Ex for each 



*The ground state doublet of mesons have light degrees of freedom with spin-parity sj l = ^ ■ 
We consider situations when only one member of this doublet contributes. It is this state that we 
denote by M. 
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FIG. 1. The integration contour C in the complex e plane. The cuts extend to Ree — * ±00. 

particle in the state X. 

While Tpvie) cannot be computed for arbitrary values of e, its integrals with appropriate 
weight functions are calculable in perturbative QCD. Consider integration of the product 
of a weight function W A (e) with T^(e) along the contour C surrounding the physical cut, 
shown in Fig. 1. Assuming W is analytic in the shaded region enclosed by this contour, we 
get 



The positivity of \{X\ J ■ a\B)\ 2 for all states X gives an upper bound on the magnitude of 
form factors mediating B decays into the ground state doublet M. 

The integral of the correlator weighted with eW&(e) eliminates the contribution from 
the ground state doublet X = M, yielding 



deW A (e) [a^T^e) a"} = ]T W A {E X - E M ) (2n) 3 5 3 (q + p x ) 



(X\J-a\B) 



(3) 



2ni Jc 



x 



2ms 




E W A (E x -E M )(27i) 3 5 3 (q + p x )(E x -E M ) 



(X\J-a\B) 



(4) 



X+M 



2vtib 
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This can be turned into an upper bound on the contribution of excited states (X ^ M) 
to the right-hand side of (0) by assuming that the contribution of multi-hadron states is 
negligible below the first excited meson state, Mi. This is true in the large N c limit, and 
experimental data available in the future on B — > D^niu, etc., decay rates can support 
(or oppose) the validity of this assumption. 
Thus, there are upper and lower bounds 



1 

2ni Jc " w 1 ^ y ' 1 AttibEm 



deW A (e) [a^T, u (e)a»] > ± (5) 



> ^- J c deW A (e) [a^T au (e)a v ] (l 



where Em x = y + \q\ 2 - ^ should be emphasized that while this upper bound (which 
yields the Bjorken bound) is essentially model independent, the lower bound (which yields 
the Voloshin bound) relies on the above assumptions about the spectrum of the final state 
hadrons X. 



Following |TJ] , we choose a set of weight functions 

A In 

that satisfy the following properties: (i) Wa is positive semidefinite along the cut so that 
every term in the sum over X on the hadron side of the sum rule is non-negative; (ii) 
Wa(0) = 1; (Hi) W^a is flat near e = 0; (iv) and falls off rapidly to zero for e > A. This 
choice of weight functions is motivated by the fact that for values of n of order unity all the 
poles of W A n ^ lie at a distance of order A away from the physical cut. As n — > 00, W A n ^ 
approaches 8(A — e) for e > 0, which corresponds to summing over all hadronic resonances 

(n) 

up to excitation energy A with equal weight. In this limit the poles of W A approach the 
cut, and the contour C is forced to lie within a distance of order A/n from the cut at e = A. 
In this case the evaluation of the contour integrals using perturbative QCD relies on local 
duality [|l^] at the scale A. In the rest of this paper whenever the weight function is not 
specified explicitly, we mean 9(A — e). 
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FIG. 2. Feynman diagrams that contribute to the order a s corrections to the sum rules. The 
heavy quark fields in the effective theory are denoted by h^ c,b ^ . The black square indicates insertion 
of the b — > c axial or vector current. 

The bounds in eq. ([|) become weaker as A is increased. However, the scale A must be 
chosen large enough that the contour integrals in eq. (0) can be performed using perturbative 
QCD, allowing the evaluation of the Wilson coefficients of the operators that occur in the 
operator product expansion for the time ordered product of currents. In practice this means 
that A must be greater than about 1 GeV. 



III. BOUNDS ON THE ISGUR-WISE FUNCTION 

The bounds stemming from eq. (|5|) are simplest to evaluate in the heavy quark effec- 
tive theory (HQET) [IT]]. One may consider the vector current in the effective theory, 
= hffl 7 M /4f), and choose a M = t> M , the four- velocity of the B. Instead of calculating 
the correlator itself, it is simpler to compute its imaginary part given by the diagrams in 
Fig. 2. In this paper we focus on the region near zero recoil, and therefore we expand the 
perturbative corrections to linear order in {w — 1). 

The ground state contribution is (M\ V 1 v^\B) = y/rriB trm (1 + w )£,( w ), where £(w) 
is the Isgur-Wise function. (Only the pseudoscalar member of the ground state doublet 
contributes to this matrix element.) In the MS scheme, using dimensional regularization 
and a finite gluon mass, m g , the inclusive expression for the correlator to order a s reads 



I, m| T( f )] = i±^ 

7T 2W 



w — 1\ , „ s 8a, , m 



SU-A +( t£ ,_i)_iin-f tf(e) 

8a, 2e 2 + m 2 r- 

+( w ~ !) q V e ~ m 9 9 ^ ~ m 9> 



+ .... (7) 



9vr e 4 

In eq. (0) the terms proportional to delta functions come from the charm quark final state, 
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and the term proportional to 6{e — m g ) arises from final states with a charm quark and a 
single gluon. The gluon mass is used to regulate an infrared divergence in the integral of 
Im [T(e)] over e, which cancels (at order a s ) between these two types of final states. 
Using eqs. (H) and (ffl) this implies the following sum rules 



1 + w 

2w 
1 + w 

2w 



1 + (w - 1 
1 



5«, 



9tt 

A(w-1) 



In. 



m Ml - m M 
3a., 



4A 2 5 
~pT ~ 3 
16a, A 



+ (w-l) 



In 



9tt A 
4A 2 5 
3 



4w 



[1 + w)'' 
Aw 



l£HI 2 + ..., 



I^HI 2 -.... 



(8) 



' 9tt \ p 

The ellipses in these equations denote positive terms whose first derivatives at w — 1 are 
also positive. The reason for positivity of the first derivative is that in the effective theory 
all excited state contributions must vanish at w = 1, and may therefore be written as (w — 1) 
times the square of some form factor. (Eq. (H) was previously obtained using a Wilson line 
approach to heavy quark interactions in Ref. |L8[]. See also fl9fl .) 

In eqs. (H), a s is evaluated at the subtraction point /i. This /i-dependence on the left- 
hand side of eq. (H) is cancelled by the /i-dependence of the Isgur-Wise function, which we 
define in the MS scheme. Differentiating with respect to w, we find the following bounds on 
the slope parameter of the Isgur-Wise function, p 2 = —d£(w) / dw\ w= i, 



1 

- + 



A 



+ 



4a, 



2 (mMi - m M) 9tt \m Ml - m M 



2A 



H In 



4A 5 



> P 2 {V) > 7 + 



4a, /5 



4 9vr V3 



- - In 



4AJ 
7?" 



(9) 



Neglecting the order a s corrections, these are precisely the Bjorken and Voloshin bounds 
discussed in the introduction. The upper bound on p 2 receives a perturbative correction of 
order a s A/Aq CD , which is very large in the A 3> A QCD limitjj] Note, however, that higher 
orders in perturbation theory do not produce additional powers of A/Aqcd- Similarly, sum 



^One is free to absorb all or part of this correction into a redefinition of A, A — > A(A), provided 
one consistently reexpresses other formulae involving A in terms of this new quantity. 
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rules involving higher moments of Im[T(e)] [f20| , |l"3f will receive perturbative strong interaction 
corrections enhanced by more powers of A/Aqcd- 

The bounds on the slope of the Isgur-Wise function in eq. (||) will have a perturbative 
series without large logarithms in its coefficients if the subtraction point \x is chosen to 
be equal to A. Note that the second term in the upper bound, proportional to A, has a 
renormalon ambiguity of order unity (since the heavy quark pole mass has a renormalon 
ambiguity of order Aqcd)- This is canceled by the ambiguity in the perturbative series in 
a s that multiplies A/(mj\/ 1 — tum), the first term of which is presented in eq. (||). 

Using a weight function W^\e) other than Wj±°\e) = 0(A — e) does not affect the lower 
(Bjorken) bound on p 2 given on the right-hand side of eq. @. Therefore, for the lower bound 
on p 2 (at order a s ), using the weight function wj±°\e) does not rely on the assumption of 
local duality at the scale A. Such corrections, however, weaken the upper (Voloshin) bound 
by adding 

4a s A ( 7r \ . , 

— T-- TT-2 , 10 

9n m Ml —m M \n sin[7r/(2n)J / 

to the left-hand side of eq. (|9|). Numerical estimates of these corrections will be given later. 

So far we have focused on the perturbative corrections to the coefficient of the lowest 
dimension operator, h^'hffl, that occurs in the operator product expansion for the time 
ordered product. Higher dimension operators are of the form (v' ■ Dph®. These yield 
corrections suppressed by powers of Aq CD /A for weight functions other than W^°\e) = 
6(A-e). 

A lower bound on p 2 including order a s perturbative QCD corrections was derived in 



Ref. It corresponds to a weight function given by the phase-space of b decay, which 

is different from those considered here. The bound in |H[ appears stronger than that in 
eq. (0) because the weight function given by the phase-space falls off faster with e. 

To zeroth order in a s and Aqcd/^c,&, the constraints in the effective theory are identical 
with bounds on the slope of the measured shape of the B — > £ v decay spectra. However, 
at first order in a s , one has to combine the above results with corrections that originate from 
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matching the full theory onto the heavy quark effective theory. [This will also eliminate the 
/^-dependence from the bounds in eq. (0).] 



IV. DIFFERENTIAL DECAY RATES 

We are interested in the form factors of semileptonic B — > Iv decays, defined as 
<DW\V»\B{v)) = h+{w) {v + vT + {v _ ^ 



y?m D m B 

ih v (w)e^e:v' a v p , (11) 



{D*{v>)\V»\B{v)) _ L ,.,.,^ a/W 



y/m D * m B 
(D*(v')\A»\B(v)) 



h Al (w) (w + l)e*»- h M (w) (e* ■v)v> M - h As (w) (e* ■ v) v'" . 



^/m D * m B 

Here V 1 = cj^b and N 1 = 07^75 b are the vector and axial currents. The kinematic variable 
w is related to q 2 via w = (m B + m 2 D( , } — q 2 ) / {2m B m D (*)) . Up to corrections suppressed by 
powers of a s {m c ^) and Aqcd l m c,bi h-(w) = h,A 2 (w) = and h + {w) = h v {w) = h^iw) = 
hA 3 (w) = £(w), where the Isgur-Wise function, £(«;), is evaluated at a subtraction point 
around m c ^. 

Experimentally the differential decay rates are measured, and are usually quoted in 
terms of the functions J r B ^ D (*)(w), defined below. At tree level, and without Aqcd/^c,6 
corrections, these functions are identical to the Isgur-Wise function, so their slopes are 
equal to that of However, at order a s additional corrections beyond those calculated 

in eq. ([]) using the effective theory arise from matching the full QCD onto the HQET. 
Corrections suppressed by powers of Aqcd/^c,;; arise from higher dimension operators in 
the HQET Lagrangian, and from higher dimension current operators in the effective theory. 

With the above definitions of the form factors, and = m D (*)/m B , the differential 
decay rates are 

— i = F f r* 3 (1 - r*f (w 2 - 1) 1/2 (w + l) 2 

aw 48 ti 6 

Aw 1 — 2wr* + r* 2 



x 



1 + 



w + 1 (1-r*) 2 



V cb \ 2 \T B ^ D 4w)\ 2 . 



i = -^—2- r 3 1 + r 2 {w 2 - 1 3/2 V/ Cb 2 ^ 

aw 48 7r d 



b^dHI 2 - (12) 
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The functions Tb^d* and Tb^d are given in terms of the form factors of the vector and 
axial currents defined in (O) as 



1 + 



Aw 1 — 2wr* + r 



*2 



w + 1 (1 — r 



*\2 



-1 



1 — 2wr* + r 



*2 



1 - r 



w — 1 

iu + 1 



+ 



W 1 n 2' 

^Ai («0 + i (tw) - ^a 3 («0 - r* h A2 (w) 



1 — r 

F B ->d{w) = h+(w) - h-(w) 



We define the "physical" slope parameters, p 2 B _, D * and p%^ D , y ia 

\Fb-+d*(w)\ = [1 - p 2 B ^ D , (w-1) + ...}, 

\Fb-+d{w)\ = \Fb^d(1)\ [1 - P 2 b^ d (w-1) + ...}. 



(13) 



(14) 



Note that £>*(!) = ^Ai(l)- Due to Luke's theorem [|] ^(1) = t/a + ^(Aqcd/to^), 
while Fb^d(X) = Vv + C(Aqcd/"Tc,6)- The quantities and 7?y relate the axial and vector 
currents in the full theory of QCD to those in HQET at zero recoil. 

The order a s corrections to the relationship between p 2 B ^ D ( t) and the Isgur-Wise function 
can be computed model independently We combine the results of the previous section with 
the order a s matching corrections |22|] taken from Ref. p3j to derive bounds on the slope 
parameters. Denoting z = m c /mb, and approximating r^*' ~ z in the order a s corrections, 
the slope of the Isgur-Wise function is related to that of J- B -*d(*) v ia 



4a 



9tt p 

Using eq. (P) this implies the bounds 



m„ a, 
+ - ' 

71 



B»i>M 27 ) + 2 d B^«.) • 



(15) 



Pb^dM > T + 



4a, 



Pb^dM < 



4 9tt 
1 



In 



a 



4A^ + V 5 



s r(a s ) 



M + 



A 



5 



(l/m) 



B - >D( * ) 2m, B - D( *' ' 



(16) 



A 



16a„ A\ 4a s , m 
1 + r- + ^ In ttj 4 



9tt A 



' 1 a -5 ( " s L„ + ^5 



9tt 4A 2 tt b -> dW 



4 2(m Ml - m M ) 

The A-independent part of the order a s corrections is contained in S B _^ D ^, 
contains the order Aqcd /m c ,i corrections to p 2 B ^ D (,)- We find 



A ~(l/m) 

2m, -B-+DW ■ 



while 5 



(l/m) 
B-^DO) 
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(ot .) _ 2(l-z)(ll + 2z + llz 2 )+24:( y 2-z + z 2 )z\iiz 
6b ^ d ' ~ 27(1 -zf ' 

(a) 2(l-2)(23 - 34^ + 232 2 ) + 12(3-3^ + 2z 2 )zlnz 

" 27(l-*)3 ' [ } 

The corrections in S B ^ M depend on the four subleading Isgur-Wise functions j| that 
parametrize all first order deviations from the infinite mass limit. These can only be esti- 



mated at present using model predictions. Using the notation of [23| we find 

= -2xi(l) + 4Xs(l) - z [2xi(l) - 4X2(1) + ^(l)] 
5 /-, \ 4 1-22 + 52 2 

6 3 3(1 — z) 

#5! = "(! + *) [2xi(l) - 4X 2 (1) + 12x' 3 (l)] + • (18) 

Here prime denotes d/dw. 

Note that the bounds in eqs. ( |16|) do not rely on the assumption that m Cj b ^> A. The 
bounds on the slope of the Isgur-Wise function in HQET (where m c ^ — > oo) hold as long as 
A is large enough for perturbative QCD to be a valid way to calculate the contour integrals 
in eq. (|5|). The values for the charm and bottom quark masses only arise in matching /)! „ w 
onto the slope of the Isgur-Wise function. 

There are several scales that occur in the bounds on p 2 B ^ D („). In the limit m c ^ ^> A 
we know how to sum the large logarithms of the ratio of scales m c fc/A that occur in the 
perturbative corrections to the bounds. First one performs the matching of p 2 B ^ D (*) onto the 
slope of the Isgur-Wise function at a scale around m c ^ (e.g., y/m c mb), then one scales the 
Isgur-Wise function down from this subtaction point to A using the anomalous dimension 
for the operator h^j^h^. Finally one applies the bound on the slope of the Isgur-Wise 
function in eq. @ at the subtraction point p = A. However, since A must be greater 
than 1 GeV, for the physical values of the charm and bottom quark masses m^/A is not 
very large. Consequently, we did not sum the leading logarithms of this ratio using this 



renormalization group procedure in eqs. (16). 
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Bjorken bound (p B _^ DM > ■ ■ ■) 
B -> D* B -> D 


Voloshin bound (p B _^ DM < • • •) 
B -> D* B -> D 


A = 1 GeV 
A = 2 GeV 


0.24 0.33 
0.18 0.27 


0.95 1.05 
1.11 1.20 



TABLE I. Upper and lower bounds on P B ^ D (*) , the slopes of the functions J r B ^ D (*) (w) at zero 
recoil, that describe the shape of the semileptonic B — > £p decay spectrum. The order a s 
corrections are included, while order Aqcd/wi C) 6 corrections are neglected. To zeroth order in a s 
the Bjorken bound is 0.25, while the Voloshin bound is 0.75. 

To evaluate the bounds in eq. (|16"D we takef] ot s = 0.3 (corresponding to a scale of about 
2 GeV), A = m Ml — m M = 0.4 GeV, m c = 1.4 GeV, and m h = 4.8 GeV. However, since we 
neglect corrections of order a s (AQCD/^ c ,fe), other values for the heavy quark masses, e.g. 
m c ~ m D = {m D + ?>m D *)/A = 1.97 GeV and m;, ~ fn B = {m B + 3m_B*)/4 = 5.31 GeV, 
would be equally valid. In Table I we show the Bjorken and Voloshin bounds at order a s 
using A = 1 GeV and A = 2 GeV, for the weight function = 6 (A — e). To obtain 

these numerical results we neglected <^_^ w . Using instead of does not affect 

the Bjorken bound, but weakens the Voloshin bound by 0.02 for A = 1 GeV, and by 0.05 
for A = 2 GeV. While the Bjorken bound only receives moderate corrections to its tree-level 
value of 0.25, the corrections to the tree-level value of the Voloshin bound, 0.75, are more 
sizable and strongly A-dependent. 

To estimate the possible size of the order Aqcd/"t. c ,6 corrections in eq. (0), we use 
the QCD sum rule predictions for the four subleading universal functions |24]|. These give 
approximately xU 1 ) = °- 3 > X2O) = -0.04, x' 3 (l) = 0.02, 77(1) = 0.6, and r/(l) = 0. Note 
that the results for depend sensitively on x'i(l)? which is only calculated to order a° 

in the framework of QCD sum rules; the other subleading form factors [x 2 ( w )> X3( w ), an d 



^Unless explicitly stated otherwise, we use these values throughout this paper. 
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rj(w)} are computed to order a s . We find that these corrections reduce p 2 B _> D * by 0.3 and 
Pb^d by 0.2. However, these values are model dependent, and the uncertainties are large. 

Using eq. (|TJ) we find that the order a s perturbative corrections predict that p\^d * s 
about 0.09 larger than p\^ D *. This prediction is affected by order Aqcd/^c^ corrections [see 
eq. flTBDl, and therefore it is not model independent. However, the QCD sum rule results 
for the subleading Isgur-Wise functions predict that p 2 B ^ D — Pb^d* * s further increased. 
Therefore, an enhancement of Pb^d compared to Pb^d* by about 0.1 — 0.2 seems quite 



likely, and a precise measurement of this difference would test the predictions of Ref. |24 
(Similar results were obtained in f25f .) 

As the above model estimates for the order Aqcd/^c,?) corrections are fairly sizable, 
one should investigate whether more reliable bounds can be derived in the full theory. The 
motivation is that in the full theory one can bound the magnitude of the physical form factors 
(equal to the Isgur-Wise function plus order Aqcd/^c^ corrections in HQET). In the next 
section we derive bounds on the {w) form factor in the full theory. For comparison, we 
give here the bounds on this form factor in the effective theory approach. The slope of 
h Al {w) at w = 1, p\ , satisfies a bound of the same form as that in eq. flT5|). The order a s 
corrections to the bounds on p\ are 



{as) _ 2(1- z)(17-Az + 17z 2 ) + 6(9-3z + 4:Z 2 )z\nz 
M ~ 27(l-z) 3 



(19) 



while the order Aqcd/^c,;; terms are 



6 (l/m) = _ 2x , (1) + ^, (1) _ z [2x , (1) _ 4xa(l) + 12x , _ 1 + Z +zr]{l) (2Q) 

Neglecting the corrections of order Aq CD /m C: b this gives 0.28 < p\ < 0.99 for A = 1 GeV 
and 0.22 < p\ x < 1.15 for A = 2 GeV. 

V. SUM RULES IN THE FULL THEORY 



In the full theory, bounds on the form factor hA 1 {w) in eq. flTT| ) can be obtained from 
eq. (HD by taking the axial current and choosing the four-vector a such that a ■ v — a ■ q — 0. 
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These bounds are expected to be less model dependent at present than those derived in 
the effective theory, as there is no uncertainty associated with the subleading Isgur-Wise 
functions (they contribute at order Aqcd/^c^ in the effective theory approach). On the 
other hand, the bounds derived in the full theory receive calculable corrections proportional 
to powers of A/m c ^ at order a s , which do not arise in the HQET approach. 

The A-dependent part of the corrections away from zero recoil are straightforward to 
compute by considering the (finite) difference between the bremsstrahlung graphs in the full 
theory and in HQET. Since we want to keep the full w dependence at order Aqcd/^c,6, we 
need to include 



l + «r 

2w , parton 



1 + W 
2w / hadron 



{w-1) A 

w 2 m r 



(21) 



Neglecting terms of order a 2 ., Aq CD /m;? b , a s (AQCD/^ c ,fe), and a s (w — l) 2 , we obtain the 
bounds 



(l + w) 2 \h Al (w)\ 2 ^ 1 
Aw 

(1 + w) 2 \h Al (w)\ 2 



w 



Aw 



> 



2w 
1 + w 
2w 



(w-1) A 



w- 



+ — 

7T 



(w-1) A 

w 2 m r 



A + (w — 1)X 
V 1 a„ 



(22) 



+ 



7T 



B + (w-l)Y 



Em! — Em. 

Here A, B, X, and Y are functions of m c , nib, and A. The term V arises from the difference 
in the start of the parton and hadron cuts, and from matrix elements of dimension-5 terms 



in the operator product expansion. It is simple to extract from ||13| . |15 
V 

where 



[w 



A A 2 3 — w Ai 3 + w , Ai + 3A 2 1 

3m h w 



w 2m c w 3 



6m c w 3 



+ 



+ ... , 



(23) 



Ai 
A 2 



1 



2m b 
1 



(B{v)\hf {iDfh^\B{v)) 



()mg (B(v)\h^ 9 -a, u G^h^\B(v)}. 



(24) 



Since E Ml — E M is of order Aqcd, all terms in V contribute at least of order AQ CD /m Cjfe to 
the lower bound on l/i^-u;)! 2 . We cannot neglect these terms, as the main motivation for 
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considering the bounds in the full theory was to eliminate the order Aqcd/^c,6 uncertainties 
related to the subleading Isgur-Wise functions in the HQET approach. 

The function A was computed in ||14|| \ and B was computed in ||15|| . Denoting d = A/m c , 
the result is 

8\ d{2 + d)[2z 2 {l + d) 2 -{3 + 2z + z 2 )} 3 + 2z-z 2 



A-- 
B = A 



2 [ i±^lnz+- ) 
1-z 3 



A 



m Ml - m M 



18 (1 + rf) 2 9 
(2 + 3d + 2d 2 ) (9 + 6,z - 3,z 2 ) + 2d 2 z 2 {8 + 7d + 2d 2 ) 
54 (1 + d) 2 
3 + 2z - z 2 ln(l + d) 



ln(l + d) 



9 d 



We find that the coefficients X and Y are 



(25) 



1 + 2, 8 

in z H — 

1-z 3 



2dz(20 + 52d + 53d 2 + 18d 3 



25 



(a.) 



d (16 + 42d + 45d 2 + 16d 3 



80 + 12z 2 



45 



45 (1 + d) 4 
ln(l + d) . 



Y = X 



AA 



A 



Al 9(l + d) 4 

d2 2 (12 + 52d + 71d 2 + AAd 3 + 10d 4 ) 
+ 45 (1 + dy 

r 16 6- lld-62d 2 -83d 3 -32(i 4 



m Ml - m M m Ml - m M I 9 18 (1 + d) A 

z(3% + 113d + 122d 2 + 41d 3 ) z 2 (34 + 119(2 + 134d 2 + 49d 3 - 16d 4 - 10d 5 ) 



45 (1 + d) 4 
15 + 38z - 17^ 2 ln(l + d) 



90(1 + d) 4 



(26) 



45 d 

In Fig. 3 we plot the upper and lower bounds on (1 + w) 2 |/i J 4i(w)| 2 /(4w) over the region 
1 < w < 1.25, using eq. (|2"2]). Over this region of w, corrections of order a s (w — l) 2 that we 
have not computed are expected to be negligible. The thin and thick solid curves correspond 
to the choices A = 1 GeV and 2 GeV, respectively. The dashed curves show the upper and 
lower bounds neglecting the order a s and Aqcd/^c,6 corrections. The dotted curves show 
the upper and lower bounds neglecting the order a s corrections, but keeping the order 
AqcdA"c,6 terms. The enhancement of the difference between the upper and lower bounds, 



Ht was first computed to order A 2 /m 2 b in [13|. 
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FIG. 3. Upper and lower bounds on (1+w) 2 \hA 1 (w)\ 2 /(Aw). The thin and thick solid curves are 
the bounds including effects of order a s and Aqcd/^c,6) anci corespond to A = 1 GeV and 2 GeV, 
respectively. The dashed curves are the bounds neglecting order a s and Aqcd/^c,6 corrections. 
The dotted curves are the bounds neglecting order a s but keeping Aqcd/"2-c,6 corrections. 

which is seen to increase with w, is dominated by the perturbative corrections. The reason 
for the somewhat larger than usual deviation of 



IMi)f 



o:. s 



7T 



1 + 2, 8\ 



(27) 



from unity in Fig. 3 is due to our choice of a s = 0.3, which gives at order a s , tja = 0.96. 

In plotting Fig. 3, we used A x = -0.2 GeV 2 and A 2 = 0.12 GeV 2 . While A 2 is well- 
determined by the B* — B mass splitting, the value of Ai is more uncertain. Changing Ai by 
±0.2 GeV 2 changes the lower bound at w = 1.25 by =f0.04. (At order Aqcd/"^c,6 the value 
of Ai does not affect the upper bound for all w, nor the lower bound at zero recoil.) We 
neglected the nonperturbative corrections of order Ag^/m 2 ^ Such corrections to the lower 
bound involve matrix elements of dimension-6 operators in the operator product expansion. 
The order AQ CD /m 2 b corrections to the upper bound [on the right-hand side of the first 
inequality in (|22l) 1 are given by |l3lJT5[ 
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FIG. 4. Upper bound on (1 + w) 2 \hj^ (w)\ 2 / (4w) including order AQ CD /m 2 6 corrections for 
different values of Ai. The curves from above to below correspond to Ai = 0, — 0.2GeV 2 , and 
— 0.4 GeV 2 , respectively. 



w 2 - 1 (Ai - \ 2 )w 2 - A 2 (3 -2w 2 ) Ai 
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2m 2 w 4 



Am 2 w 5 



X 2 w 2 Ai + 3A 2 ( 2 2 . o x 
+ — — n — - [ow z + 2z) 



12m 2 w 3 



(28) 



For the central values of A and Ai used throughout this paper, these corrections affect the 
upper bounds by only about —(0.02 — 0.03). The size of this correction is sensitive to the 
value of Ai, while it is largely independent of A. In Fig. 4 we plot the upper (Bjorken) bound 
on l/i^w)! 2 including the order AQ CD /m 2 fc corrections in eq. ( |2"8"D for A = 1 GeV and three 
different values of Ai, Ai = 0, -0.2 GeV 2 , and -0.4 GeV 2 . Changing A by ±0.1 GeV from 
the value we used (i.e., 0.4 GeV) affects the curves plotted in Fig. 4 by about =f0.01. If 
Ai is known more accurately in the future, this upper bound may constrain h^l) to be 
somewhat below t\a- 

A larger value of A increases our confidence in the validity of using perturbative QCD to 
evaluate the time ordered product of weak currents. In going from A = 2 GeV to A = 3 GeV 
the upper bound (the upper thick solid curve in Fig. 3) is increased by about 0.01 fairly 
independently of w over the region 1 < w < 1.25, while the lower bound (the lower thick 
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solid curve in Fig. 3) is decreased by 0.06 at w = 1 and by 0.07 at w = 1.25. Clearly the 
lower bound is considerably more sensitive to the choice of A than the upper bound. 

Converting the bounds in eq. ( p2"D into constraints on the slope parameter of the h Al form 
factor at zero recoil, p\ ± , is not straightforward. The upper and lower bounds on l/i^w)) 2 
do not meet at zero recoil, and therefore a bound on p Ai can only be derived with some 
smoothness assumption. We bound p\ by assuming that h Al (w) is linear over the region 
1 < w < 1.25. Then the plots in Fig. 3 (which neglect terms of order AQ CD /m 2 fe ) imply the 
bounds 0.44 + 4[h Al (l) - r] A ] < p\ < 1.19 + 4[/u x (l) - rj A ] and 0.39 + 4{h Al (l) - rj A ] < 
p\ < 1.36 + 4[/i J 4 1 (l) — t)a\ for A = 1 GeV and 2 GeV, respectively. Recall that h Al (1) — f] A 
is of order AQ CD /m 2 fc . The increase in the lower bound compared to the 1/4 at zeroth order 
is mostly due to the terms proportional to A/m c in eqs. (|22|), and hence it is sensitive to 
the value of A we choose. It is interesting that this correction has the opposite sign than 
the QCD sum rule results in the effective theory, which predicted that order Aq C d /^ c ,6 
corrections lower the values of the bounds. 

The smoothness assumption used to derive bounds on the slope parameter p 2 Ai from 
Fig. 3 can be justified model independently using the parameterization of Ref. [F2(J. This 
work also lets us extend the bounds at small w presented in Fig. 3 to larger values of w. The 
A = 2 GeV bounds imply upper and lower bounds at w = 1.5 of 0.86 > h Al (1.5) > 0.36. 

The bounds presented in this section can also be used for unpolarized A& — > A c £ v decay 
if certain replacements are made. The form factors for the matrix element of the axial 
current, G i: are defined by 

(A c (v')\A"\A b (v)) = u(v') [G^ + G 2 v» + G 3 v'»] 7s u(v) . (29) 



Bounds on G± are obtained by replacing the left-hand side of eqs. ( p2|) by G\{w + l)/{2w). 
The quantities A, Ai, and A2 that appear on the right-hand side must now be interpreted as 
arising from A^ matrix elements. They are simply related to the corresponding quantities 
in the B meson case P7] , A2(A{,) = 0, 
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Ai(A 6 ) = \i{B) + 2m c [(m Ab - m B ) - (m Ac - m D )]/ (1 - z) , 
A(A 6 ) = A(S) + m Ab -m B + [Ai(A 6 ) - Ai(£)]/(2m b ) . 



(30) 



VI. ORDER a 2 /3 CORRECTIONS AT ZERO RECOIL 

At zero recoil we can rewrite the upper and lower bounds on the Ua x form factor in 
eq. as 

a s (A) A , o2(A) 

r 2 



|/l Al (l)| 2 < 7^ + A! + ^ , (31) 

7T 7T Z 



IWUI' >r,l+^>A 1 + ^>P A 2 



«.(A) D , o2(A) 



77 7T 2 



These bounds also receive nonperturbative corrections of order AQ CD /m 2 b . Such corrections 
to the upper bound are given by 

-^ + ^^( , l + ^ + * 2 V 02) 



m 2 4m 2 ^ 3 

The order AQ CD /m 2 6 corrections to the lower bound on 1/^(1)1 2 depend on additional 
quantities (besides Ai and A2) that parametrize matrix elements of dimension-6 operators 
in the operator product expansion. 

The corrections in tja in eq. (|31|) arise, at the parton level, from the final state X = c. 
Except for an infrared renormalon ambiguity, tja only depends on physics associated with 
the scales m Cjb . It has been calculated to order a 2 /3 [p8| . Explicitly, 



a s (y/m b m c ) / l + z 8\ _ a 2 s (^/m b m c ) 5_ S l + z 44 

vr U-z 112 3/ 7T 2 ft 24U-z ^ 15 



The full order a 2 expressions for t]a and 77^ are also known approximately p9| . For t)a the 
order a 2 @ terms dominate the a 2 correction, while for rjv they do not. 

The corrections in Aj and Bi originate from final states X that contain a charm quark 
plus additional partons, e.g., eg, cqq, etc. These corrections are suppressed by powers of 
A/m Cjb . For convenience, we evaluated the arguments of the strong couplings in eqs. Q5TJ ) 
and (|33"D at two different scales ^m c m b for the series in t)a and A for the other terms). 
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FIG. 5. The functions A\, B\, A 2 , and B 2 denned in eq. (|3l| ) as a function of A. The dotted 
curve is A\, the dash-dotted curve is B\, the dashed curve is A2, and the solid curve is B^- 

Of course, it is possible to evaluate both series at the same scale using the QCD /3-function 
pOfl . The functions A\ and B\ are given in Sec. V, A 2 was computed in Ref. |14|, our result 
for B 2 is given below. These order a 2 (3q corrections are relatively simple to compute due 
to a relation between the rif dependent part of the order a 2 terms and the order ot s result 
with a finite gluon mass [|3T]. The calculation is simplest in the so-called ^-scheme, but we 
present the results in the usual MS scheme. To leading order in d = A/m c 

d 2 ( 2 „\ „ 1/13 \ „ d 2 ( 4 9 \ 

2 „ 2 t d 2 ( 2 9 \ 
B 1 = - Ax , B 2 = -A 2 1 + - z + z 2 . 

3 ' 3 54 V 3 J 

To all orders in d, A x and B\ are given in eq. (p5[) , and A 2 and B 2 are computed numerically. 
In Fig. 5 we plot the functions A%, B\, A 2 , and B 2 as a function of A. The dotted curve is 
Ax, the dash-dotted curve is Bi, the dashed curve is A 2 , and the solid curve is B 2 . 

On the right-hand sides of eqs. ( j3~H ) the renormalon ambiguity of order AQ CD /m^ fe in rf A 
is cancelled by that in the series A\. Therefore, in considering the size of the a 2 corrections 
relative to those of order a s , it seems reasonable for the upper bound to take all terms on the 
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(34) 



right-hand side of eqs. ( [31]) together |32 |. Although A 2 is approximately as big as A\ (over 
the region A < 2GeV), since the order a 2 /3 term in rj A partially cancels against A 2 , the 
perturbative expansion for the upper bound seems reasonably well-behaved.Q For example, 
with Aqcd = 200 MeV and A = 1 GeV, one finds neglecting terms of order AQ CD /m^ fc that 

\h Al (l)\ 2 < 1 - 0.073 - 0.019 
+ 0.013 + 0.017 
= 1 - 0.060 - 0.002. (35) 

The last two terms in the first line are the order a s and a 2 (5q corrections to 7] 2 A) while the 
second line contains the terms proportional to Ai and A 2 . For values of A between 1 GeV 
and 5 GeV the cancellation of order a 2 /3 terms persists. 

For the lower bound, the order a 2 (3$ term in the perturbative expansion of the term 
proportional to A/(rriM 1 — ^m), which originates from the first moment of the time ordered 
product of weak currents, is about as large as the order a s term over the region A < 5 GeV. 
With A QCD = 200 MeV and A = 1 GeV, the lower bound in eq. @ is 

\h Al (l)\ 2 > 1 - 0.060 - 0.002 
- 0.019 - 0.022 
= 1 - 0.079 - 0.024. (36) 

The terms in the first line arise from rf A and the series A iy while the second line contains 
the terms proportional to B\ and B 2 . Note that the coefficient of the Bi terms depend on 
the mass of the first excited state. 



**Ref. was less certain that the upper bound on |^i J 4 1 (l)| 2 has a well-behaved perturbative 
expansion. This is mostly due to the fact that in [l4| the behavior of the perturbative series arising 
from the terms Ai (which contain all the A-dependence) was discussed without combining it with 
that from r] A . 
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VII. CONCLUSIONS 



In this paper we studied corrections to the Bjorken and Voloshin sum rules on form 
factors of semileptonic B — > tv decays. In the heavy quark effective theory we derived 
upper and lower bounds on the Isgur-Wise function, and on its slope at zero recoil. Matching 
the full theory onto HQET, we translated the bounds in the effective theory into bounds 
on J r B _ tD (*), the shape of the measured B — > spectrum. The results in Table I show 
that while the corrections to the Bjorken bound (lower bound on ^ „ w ) are small, the 
corrections to the Voloshin bound (upper bound on p 2 B ^ DM ) are sizable. The reason is that 
perturbative corrections to the Voloshin bound are enhanced by A/A. Therefore, even if 
experimental data would settle around p B ->D(*) slightly above unity, that would still not be 
a problem for the theory to accommodate. 

The bounds derived in Sec. IV are affected at order Aqcd/^ c ,6 by corrections that are 
parametrized by four subleading universal functions, and are not known at present model 
independently. Therefore, we also studied the sum rule constraints on the Ha x {w) form 
factor in the full theory. In this approach, all Aqcd/^c6 corrections to the sum rules are 
parametrized by the three matrix elements, A, Ai, and A2. Bounds on the Hai form factor 
are shown in Fig. 3, and with some assumptions we found that its slope parameter satisfies 
0.4 < p\ < 1.3. 

These bounds on h^ x in the full theory of QCD can also be related to bounds on p 2 B ^ DM ■ 
Using eqs. Q17D and ( |19|) we find that the order a s corrections imply p 2 B ^ D * = p\ —0.04 and 
p 2 B ^ D = p\ r +0.05. These relations receive order Aqcd/^i c ,6 corrections. Due to heavy quark 
spin symmetry the order Aq CD /m^ b difference between p 2 B ^ D (*) and p\ is independent of 
the subleading Isgur-Wise function \i- Then eqs. ( pT8|) and (|20|) , together with the QCD 
sum rule predictions for the subleading Isgur-Wise functions, imply that order Aqcd/"t. c ,6 
terms reduce p 2 B ^ D * — p\ x by 0.11 and reduce p 2 B ^ D — p\ by 0.01. (The uncertainty in these 
predictions will be reduced if the form factor ratio R2 can be measured precisely.) 

One of the largest uncertainties in the sum rule predictions (especially in the Voloshin 
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bound) is related to the numerical values of A and mj/j —tum- In this paper we used mm, — 
rriM = A = 0.4 GeV, motivated by the experimentally measured Di — D* mass difference, 
and by the extraction of A in p3]| . However, the uncertainties in this determination of A 



are sizable, and precise experimental data on other inclusive processes is needed to extract 
the value of this quantity more reliably [[53J]. Taking iriMx — n^M — ^d x — m D* may also be 
misleading. Besides the possibility of sizable decay rates into non-resonant final states |35| 



(discussed in Sec. II), there is probably a doublet lighter than the {Di^D^} heavy quark 
spin symmetry doublet, that is of order 100 MeV or more broad. The spin-parity of the 



light degrees of freedom in the {Di^D^} doublet is sj l = | + , and so the D\ is 1 + , while 
the D\ is 2 + . However, light degrees of freedom with sj l = \ + yield a doublet of + and 
1 + states. These can decay into D^'tt in an s-wave, and so they should be much broader 
than the {.Di,/^} that can only decay in a <i-wave. (An s-wave decay amplitude for the 
D\ is allowed by angular momentum conservation, but it is forbidden by heavy quark spin 
symmetry |5B| .) A reduction in — itlm would further weaken the Voloshin bound. 

Some improvements in this paper are possible. We focused on the region near zero recoil 
because it is important for the extraction of |V^,|. It would be straightforward to calculate 
the terms in eq. ( F2"2"D of order a s (w — l) 2 . Then Fig. 3 could be extended over the full 
kinematic range 1 < w < 1.5. Some uncertainty in the sum rules arises from the order 
a 2 s corrections. The part of these corrections proportional to the one-loop /3-function have 
been computed at zero recoil. For the part of the lower bound involving the first moment 
of the time ordered product of weak currents these corrections are as big as the order a s 
corrections, unless A is quite large. It should be possible to compute the a 2 s (3q corrections 
away from zero recoil, and also to the bounds on p 2 B ^ D ^) considered in Section III. 
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